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Abstract 

This work develops an e-uniform finite element method for singularly perturbed 
boundary value problems. A surprising and remarkable observation is illustrated: By 
moving one node arbitrarily in between its adjacent nodes, the new finite element 
solution always intersect with original one at fixed point. Using this fact, an effective 
e-uniform approximation out of boundary is proposed by adding one point only in the 
grid that contains boundary layer. The thickness of boundary layer is not necessary to 
be known from priori estimation. Numerical results are carried out and compared to 
Shishkin mesh for demonstration purpose. 
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1 Introduction 



This paper is concerned with linear Galerkin finite element method for singularly perturbed 
boundary value problems (BVPs). Consider an one-dimensional BVP problem 

-eu" - bu' + cu = f, u(0) = u(l) = 0, re e [0,1]. (1) 

For simplicity, let b < 0, c > 0, and < e <C 1 are constant such that not both b and c are 
0. If b > 0, by using substitution w(x) = u(l — x), it reduces to the case with b < 0. If 
b = 0, c > 0, equation ((TJ) is said to be a reaction diffusion equation. If b < 0, c = 0, equation 
(JTJ) is the so-called convection diffusion equation. All the results presented in this paper can 
be readily generalized to smooth and non- vanishing functions of b(x) and c(x). 

If the exact solution u(-) of (0) is "bad" in the sense that ||w"||oo is n °t bounded uniformly 
in e, the standard finite element method (FEM) generates huge errors through the whole 
domain. Typically, it is caused by a small interval of width 0(e) (called boundary layer), in 
which u" rapidly changes. 

To overcome the difficulties in the singular perturbation, it is desirable to put more grid 
points near the boundary layer or stablize the appoximation methods. Streamline diffusion 
finite element methods (SDFEM), upwinding FEM, Bakhalov grid, Shishkin grid, and many 
other such schemes are extensively studied in the context of singularly purturbed problems 
since 1970s, see [TTJ IH El IJ3 Hj- Among them, Shishkin grid became popular due to 
its simple structure and high accuracy. The Shishkin mesh was first introduced in finite 
difference methods and has been discussed in jS]; the reader is referred to a survey article 
[13 j for further details. A typical Shishkin mesh is to construct n + n grid, which is indeed n 
uniform grids in boundary layer plus n uniform grids out of boundary layer. By this method, 
the approximation provides e-uniform accuracy. But they require a priori estimation in order 
to determine the thickness of the bounded layer. On the other hand, it makes error analysis 
more complicated, since the errors from boundary layer affect the solution in the entire 
domain. Therefore, if an approximation can be stabilized and e-uniform by simply adding 
one point to original n grid, it deserves to be worked out. 

In this work, we focus on FEM solutions of (0) by starting with an interesting observation. 
Given a grid T n = {0 = Xq < ■ ■ ■ < x n < x n+ i = 1}, we add m points arbitrarily in 
[x n , 1], denoted by T n+m . Then FEM solutions on T n and T n+m intersect each other at a 
fixed point in each interval of out of boundary layer, that is, the locations of intersection 
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in each interval is indpendent of m and the distribution of added points, see Figure ^-a, 
Figure OJ-a, and Figure EJ This directly implies that the accuracy on those intersections 
are as good as FEM solutions on the grid T n+m by choosing m — > oo, denoted by 7" n +°°. 
Provided that the boundary layer is covered by [x n , 1], the above observation gives the start 
point of e-uniform approximation. In lieu of interpolating these intersections, we present a 
better way to obtain an e-uniform approximation. By adding one point s\ G (x n , 1) with 
\si — x n\ — 0(e) or \§i — x n \ = 0(y/e), the interval [x n , Si] block the error impact from 
boundary layer completely. The theoretical result shows that the FEM solutions with grid 
T n U {§i} in [0, x n ] is the same as the FEM solution of 

— ew" — bw' + cw = f, w(0) = 0, w(s\) = u(§i), x G [0, s\], (2) 

where u(s\) is exact solution at the point of si, and w"(-) is uniformly bounded. This enables 
us to use all kinds of standard FEM error analysis in [0, s\], no matter how huge errors are 
generated in [s"i, 1]. Therefore, the FEM errors in [0, x n ] has the accuracy of 7"™+ 1 + 00 ) which 
is clearly better than Shishkin mesh T n+n in both accuracy and computing cost. Another 
advantage is: One need not know the thickness of boundary layer, since s\ is not necessarily 
in boundary layer. 

The rest of this paper is arranged as follows. Section 2 begins with the model and 
notation. Section 3 proceeds with the observation on intersections of a family of FEM 
solutions. Section 4 presents an e-uniform FEM, which can isolate the boundary layer. Some 
auxiliary results are included in Section 5. Section 6 displays some numerical experiment 
results, including solutions of convection-diffusion equation, reaction-diffusion equation, and 
Green function. Finally, we close this paper with further remarks. 

2 Formulation 

Let H 1 = {v,v' G L 2 }, and Hi = {v\v G #\t;(0) = v(l) = 0}. The weak solution of © is 
a function u G Hq, satisfying 

a(u,v) = (f,v),VveH 1 , (3) 

where (•, •) is L 2 inner product, and a(u, v) = e(u',v') + b(u,v') + c(u,v). 
For a positive integer n > 2, let T n be an arbitrary grid of the form 

T n = {xi\0 = x < x 1 < ■ ■■ < x n+1 = 1}, (4) 
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and let hi = — By 4>i{ x )i we denote the nodal basis function at Xj for 1 < i < N by 

( X — Xi-\ 



h; 



if x G [x^Xi] 



if x G [xi, x i+ i] (5) 



otherwise. 



The finite element space is defined by V n = {v n \v n = Y17=i ,y ? < M ;r )}- The finite element 
discretization of (jHJ) is to find u n G V n such that 

a(u n ,v n ) = {f,v n ),Vv n G V n nHl (6) 

Existence and uniqueness of u n can be found in [2] and references therein. Now we denote 

n 

u" = (7) 
i=i 

Rewrite (0) as 

^ fy) = (f, <f)j),j = 1, 2, . . . n. (8) 

2=1 

Let i be an n x n matrix with 

dij = a(4>j , 4>i) . (9) 
Detailed calculation leads to further specific form of 

/ 1 1 N C /7 j \ 

=e(— + T — ) + + 

e b c 

QiA-i — —-, — I" 77 + -^hi 
hi 2 6 

£ & c 

tt* 7-1-1 = ^ I — 7^7 

' 2 6 

a ?ii = 0, if \i -j\>2 



(10) 



Let U n = «,...,<)' and F = ((/, 0i), ...,(/, </>„))' be column vectors. Then, © is 
equivalent to the linear system of equations 

AfT = F. (11) 

Typically a FEM solution of a singularly perturbed BVP problem has boundary layer in 
a small interval (associated with e) of rapid variations of u". Throughout this paper, unless 
it's explicitly mentioned, we assume solution u of (JTJ) has a boundary layer at x — 1 and x n 



4 



is located outside the boundary layer. This is reasonable assumption due to the very short 
interval of boundary layer depending on < e <C 1. All the results below can be obtained 
analogously for any boundary layer located in [0, 1]. 

Let T n+m = T n U {si, . . . , s m }, where x n < s 1 < ■ ■ ■ < s m < x n+ %. Denote the 
nodal basis functions on T n+m by . . . , n _i, 4> n , 4> S1 , . . . , 4> Sm }, where n and (p Si are 
nodal basis for x n and Sj, respectively. Note that the first n — 1 nodal basis functions 
of T n+m are exactly the same as those of T n . Let V n+m be the function space with ba- 
sis {0!, . . . , n _i, n , (j) sx , . . . , <p Sm }. It is obvious that V n C V n+m . Write u n+m , the FEM 
solution of (JTJ in V n+m , as 

n— 1 m 

u n +m = uftonfr + u n n +m ^ n + < m 4> Si . (12) 

i=l i=l 

In the next section, we fix T n , and start with observation on the intersections of u n and 
u n+m £ Qr diff eren t q-n+m^ p Qr convenience, we use Qi G u n r\u n+m to denote the intersetion of 
u n and u n+m in the interval Xj), and by and we denote x- and y- coordinate 

of Qi respectively. The result shows that the intersections {Qi : 2 < i < n} are independent 
of m and distribution of Sj. Therefore, by adding only one point {si}, we can compute 
{Qi 6u"n M n+1 }, and the accuracy of Qi has the same accuracy as -u n+ °°. 

3 Intersections of u n and u n+m 

Theorem 3.1. Fix T n . By adding one point Si G (x n , 1) arbitrarily, we obtain new grid 
T ra+1 . Then the intersection Qi of u n and w n+1 in the interval is independent of 

the choice of S\ for any i = 2,3, ... ,n. That is, those coordinates of intersections do not 
depend on the choice of Si G (x n ,x n+ i). 

Proof. Analogous to (jHJ), we have a system of linear equations with respect to {u™ +1 ,i = 
l,...,n; m^ 1 }, given by 

n-1 

^ < +1 a(0 4 , + < +1 a(0 n , 0,) + < +1 a(0 Sl , 0,) = (/, 0,), j = 1, 2, . . . , n - 1, (13) 
i=i 

n-1 

J]< +1 a(0 4 , 4> n ) + < +1 a(0„, n ) + < +1 a(0 Sl , n ) = (/, n ), (14) 
i=i 
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and 

n-l 

5>7 +1 a(&,& x ) +< +1 a(0 n ,0 Sl ) + < +1 a(0 Sl) Sl ) = (15) 

8=1 

Note that for 1 < j < n — 1, a(0 n , 0-,) = a(0 n , 0j) and 0(0^, 3 -) = 0, and (fT3|) leads to 

^< +1 a(0;,0,) = (/,0,), j = 1,2,..., n-l. (16) 

i=l 

On the other hand, for 1 < i < n — 1, a(0j, n ) = a(0j, n ), and (fT^j) yields 

n-l 

^< +1 a(0,,0 n ) +< +1 a(0 n ,0 n ) = (/,0 n ) -< +1 a(0 Sl ,0„). (17) 

i=l 

For 1 < i < n — 1, a(0i, S1 ) = 0, so it follows from (fTo'j) 

< +1 a(0 n ,0 Sl ) = (/,0 S1 ) -< +1 a(0 Sl ,0 Sl ). (18) 

Let p = (1 — s)/h n+ i. Observe <p n = (fi n + pcp si . Combining two equations above according 
to (fT7 )l +E>* ([T3|l . we have 

n-l 

^< +1 a(0 i; n ) + < +1 a(0„,0n) = (/,0 n ) -< +1 a(0 Sl ,0 n ). (19) 
j=i 

Hence, 

n 

^< +1 a(0 i? n ) = (JM -< +1 a(0 sl ,0 n ) +p< +1 a(0 Sl ,0 n )- (20) 

i=l 

Let £/ n+1 = . . . , < +1 )' be a column vector with length n. By (fTCjl and (f27Ijl . 

Af/ n+l = ^ ^ 

where F is a column vector with left-hand side of (jlfij) and (f2T)|) as elements. Subtracting 
(J2H) from (HH), 

A([T - U n+1 ) = F-F (22) 

Notice that F — F = C s ^e n) where e n = (0, . . . , 0, 1)' is a vector with length n, and C s> i = 
w" 1 +1 a(0 si , 4> n ) — pu™ +1 a(0 si , 4> n ). Note that C s> i is a scalar depending only on si, since 
term in C Si i is completely determined by Si. Therefore, 

U n - f/ n+1 = C s>1 A- l e n (23) 

The last equation tells us every — increases or decreases by the factor C Si i uniformly 
in i. Using elementary similar triangle properties, we prove the result. □ 
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Remark 3.2. If uf — w" + and — w™^ have opposite sign, then u n and u n+l have 
intersection in (xj,£j +1 ). Notice that A~ 1 e n in (|23|) is FEM solution of green function of 
operator A. It is very common that FEM solution of green function intersects x-axis in 
each grid. Intuitively, this explains why u n and u n+l intersect each other in every grid in 
most cases. Later we will present the criteria to be used for identifying the existence of 
intersections, see Lemma f5. II Moreover, if there is no intersection in some interval (xi,Xi + i) 
for a choice of s\, then there will be no intersection for any choice of s\. 

Theorem 3.3. Fix T n . Let T n+m = T n U {s x < s 2 < ■ ■ ■ < s m }, where Sj G (x„, 1). 
Then the intersection Qi of u n and u n+m in the interval (xi-i,Xi) is independent of m and 
distribution of {s^} for any i — 2, 3, . . . , n. 

Proof. Let V n+m be a function space with nodal basis functions {0i, . . . , 0„_i, <fi n , <f) si , . . . , 4> Sm \ 
on T n+m . Analogous to (fTBj) . we have 

n 

J]< +m a(0 l ,^) = (/,^), j = l,2,...,n-l. (24) 

Since V n+m D V™, there exists a linear combination n = n +X]^li P«0si f° r some pi,p 2 , ■ ■ ■ ,p m G 
[0, 1]. Applying similar arguments as that of Theorem 13.11 we obtain 

n m 

J]< +m a(^,0 n ) = (f,<p n ) + J2(P< +m ~ u^ m )<4>s v 4>n)- (25) 

i=l i=l 

Define C Sjjra = ^™ =1 (w" i +m — PiW^ +m )a(0 Si , n ). Using exactly the same argument in (|2"3*j) . we 
have 

f/ n - U n+m = C s , m A- x e n . (26) 

Hence, the result follows. □ 

Corollary 3.4. Fix T n . Let T n+m = T n U {s x < s 2 < ■ ■ ■ < s m }, where s< G (0,xi). 
Then the intersection of w n and M n+m in the interval (xi-i,Xi) is independent of m and 
distribution of {sj} fixed for any i = 2, . . . , n. 

Proof. We rearrange the order of the index from {0, 1, 2, . . . , n, n + 1} to {n + 1, n, . . . , 1, 0}, 
and change the coordinate linearly from [0, 1] into [1,0]. Using the same line of argument as 
that of Theorem 13.31 the result holds. □ 
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Corollary 3.5. Fix T n . Let T n+m = T n U {s 1 < s 2 < ■ ■ ■ < s m }, where s< G (x k ^,x k ) 
for some 2 < k < n. Then the intersection of u n and u n+m in the interval is 
independent of m and distribution of {sj} for any i G {2, . . . , n} \ {k}. 

Proof. This is straight forward result from Theorem 13.31 and Corollary 13 . 41 □ 

4 An e-uniform Approximation u n+l in [0, x n ] 

In the previous section, by arbitrarily choosing a point S\ G (x n , 1), we can determine 
Qt G u n fl in each interval, and the result shows Qi has the same accuracy as that 
of M n+ °°. In this section, by choosing appropriate s± G (x n ,l), we obtain u n+1 , which has 
e- uniform accuracy in [0, x n ]. This will automatically imply that Qi has e- uniform accuracy, 
since Qi G u n+1 . For simplicity, we slightly abuse notation: Let a n ^ = a((p s ~. , (p n ) without 
confusing. 

Lemma 4.1. There exists Si G (x n , 1), such that, a n5l = for T n+1 = {x < x\ < ■ ■ ■ < 

X n < S\ < X n+ \\. 

Proof. By f)l()jl . to establish the desired result, it is equivalent to prove that there exists 
< h§ 1 < 1 — x n , satisfies 

-i-i + f^.O, (27) 
where h Sl = §i — x n . By eliminating the denominators in the equation (|2*7jl . we have 

ch\ - 3bh Sl -6e = 0. (28) 
— 2e 

If c = 0, then b < 0, and hg, = — — > 0. If c ^ 0, then the determinant of is 

o 

96 2 + 24ec > 0. Write h§ 1 using quadratic formula, 



36 + + 24ec /6e , N 

< h„ = < y/-. (29) 

Thus, h tl = 0{e) if c = 0, and h h = 0{^e) if c ^ 0. □ 

Remark 4.2. The essence of Lemma [4. II is to find such a hg 1 with a^^^ =0. If b and c are 
not constant, we can compute the formula for hg 1 involved with integrals. It is also possible 
to find it by discretizations. 
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Theorem 4.3. Given T n , take T n+1 and §i as in Lemma HHJ Use u n+l to denote the FEM 
solution on T n+1 of ([Tfl. Consider another BVP problem 

-ew" -bw' + cw = f, w{0) = 0,w(si) = u(si), (30) 

where «(■) is solution of Use w n to denote the FEM solution of (fHUJ) on T n+1 \ {1}, then 

w(x) = u(x), Wx G [0,si], (31) 

and 

u n+1 (x) = w n {x), Vx G [0,x n \. (32) 

Proof. Note that , W2 +1 > • • • > ^n +1 ) * s a solution of the system of linear equations 
ai,i-i<+i + a M < +1 + a iA+1 uj£ = (/,&) i = 1, 2, . . . , n - 1 

Let w n = Ya=i Wi&i+Wn&n+Wg^, where 0r =0^| [o5i] . Then (u>™ for i G {1, 2, . . . , n, si}) 
is a solution of the system of linear equations 

a^-iiuJL.! + + a M+1 w™ +1 = (f,(pi) i = 1, 2, . . . , n - 1 

a^n-xiy^ + a n , n u£ = (/, n ) - cin^w? (34) 

w\ = u(h)- 

The solutions of (|3^j) and (|3*4*jl are precisely the same, since a n ^ x = 0. □ 

Remark 4.4. From Theorem 14. 3| we can separate the boundary layer by adding point s\ G 
(x n , 1). Therefore, it is equivalent to solve non- singularly perturbed BVP problems by the 
FEM, and all general FEM error analysis works well without effected by boundary layer. For 
example, if {0 < xi < ■ ■ ■ < x n } is uniform mesh in [0, x n ], then ||w n+1 ||oo,[o,:r„] is bounded 
by IKIIoo^o^+CHVi)]^ 2 , and IK / |L,[o,*„+o(^)] is e-uniformly bounded. On the other hand, 
add m points in (Si, 1), denoted by 7" n+1+m . Use 'j~ n+1+00 to denote the grid which is almost 
dense in [s%, 1]. Use to denote the FEM solution of JTJ) on f n+1+ °°. Then, u n+1 is 

exactly the same with £ n + 1 +°° on [0,x ra ]. 

5 Auxiliary results 

Recall A is an n x n matrix with 4>i), and e n is (0, . . . , 0, 1)' of length n. Let 

be matrix replacing zth column of A with e n . 
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Lemma 5.1. Fix T n . Let T n+m = T n U {s 1 < s 2 < ■ ■ ■ < s m }, where Sj e (x n ,x n+1 ). Then 
w n and u n+m have their intersection Q« in the interval for some 2 < i < n if and 

only if 

det[44i-i] < 0, (35) 

and the coordinates of Qi is given by 

Qi = ( — 7-7^ + —77^-1, —7T U ? + —TT u ti ) > (36) 

\n + 1 Ti + 1 Ti + 1 7-j + 1 / 

where r { = \ detA^/ det Ai\. 

Proof. To obtain Qi, we apply T n and T n+m to Theorem 13.31 Using Crammer's rule in 
(l2*oTl . we obtain 

Hpt 4 • 

<_<+m = a ' m S' i = l,2,...,n, (37) 

Therefore 

det A 2 _ 2 n (3g) 
uU-utr "det A_ x ' ^" 2 '---' n - ^ 

A necessary and sufficient condition to have an intersection is (uf — u'i +m ) / (m™+i — M™+i m ) < 0. 

This proves 1)35)) . Using similar triangles, (j36| follows. □ 

It is very common to have oscillation in finite element solution, and we can use Lemma l5~Tl 

to verify its behavior, see Remark 13.21 The following theorem is a direct consequence of using 

Shishkin mesh. 

Theorem 5.2. Assume T n is a uniform grid in [0, 1] satisfies condition (|35jl . and the bound- 
ary layer is at x = 1. Then 

max \u n (x(Qi)) -u](x(Qi))\ <Cn~ 2 } (39) 

l<i<n— 1 

where C is independent of e. 

Proof. We put m = 0(n) grid in (x n , 1), so that T n+m forms Bakhvalov grid or Shishkin 
grid. The uniform convergence of u n+m on T n+m is well known (see [3 El E3) as 

||w n+m - < +m ||oo < CrT 2 (40) 
Also, we have Qi G u n+m PI «" by Corollary 13.51 So 

K +m (x(Qi)) - < +m (x(g,))l < Cn- 2 . (41) 
Note that w" +m |(o,:r n ) = u l\{o,x„)- Thus, the theorem holds. □ 
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Remark 5.3. From the result of Theorem 15 .2\ we have estimation of 0(h 2 ). In non-uniform 
case, we can obtain an error bound 0(h) directly from j^j. 

6 Numerical Results 

In this section, we present several examples. The first is a convection diffusion equation, the 
second is a reaction diffusion equation, and the last one is a Green function. 

Example 6.1. Consider the convection-diffusion equation: 

-eu" + u' = x, u(0) = u(l) = 0. (42) 

The exact solution is 

U = *(l + S )-(2 +S )( 1-e-V. )■ < 43 > 

The solution u(-) has a boundary layer at x — 1, and is nearly quadratic outside the boundary 
layer. 

First, we use the linear finite element method on two different grid T 15 and T 15+1 for 
e = 10~ 3 , where T 15 is a uniform mesh on [0, 1] with 16 intervals, and T 15+1 is a modified 
T 15 with one point added at the center of the last interval. The intersections of finite element 
solution u 15 and m 15+1 are almost on the interpolation of exact solution shown in 

Figure [TJ-a. 

Second, we use the grid T 15+1 to compute for the same e, where T 15+1 is modified 
from T 15 by adding one specific point s\ G (x n , 1) with Si — x n = 2s, see Lemma f4.ll The 
finite element solution u 15+1 is almost overlapped with interpolation of interpolation of exact 
solution u} 5 in [0,x n ], as seen from Figure ^-b. This verifies Theorem 14.31 

To compare with the well-known Shishkin mesh, we construct T s n+n , which divides both 
[0, 1 — 9] and [1 — 6, 1] into n equidistant subintervals, where 9 — {|, 2e ^ 2n }. u™ +n is used to 
denote the FEM solution on T™ +n . Tabled shows the maximum norm of u^ +l — u n+1 and 
— u™ +n in [0, x n ]. Apparently, both u n+1 and u" +n has e-uniform accuracy. However, 
yn+i j-jg^gj. accuracy than u™ +n by using less grids. The reason is that u n+1 is completely 
isolated from the impact of errors from boundary layer; see Table ^ This also verifies 
Theorem 14.31 
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e = lo- 



rn 



-10 



\n n+l — i"i n + l \\ 

\Uj —U ||oo,[0,x„] 



\ii n+1 — ii n + n \\ 

\Uj — U s ||oo,[0,a: n ] 



U/ n+1 _ f,n+l II r , |L, rl + 1 _ n,n+n || 

Pj « ||Do,[0,x n ] \\ u l ~ u s l|oo,[0,x„ 



4 


6.663e-003 


1.117e-002 


6.667e-003 


1.117e-002 


8 


2.054e-003 


1.567e-003 


2.058e-003 


1.569e-003 


16 


5.734e-004 


3.480e-004 


5.767e-004 


3.500e-004 


32 


1.498e-004 


8.384e-005 


1.530e-004 


8.569e-005 


64 


3.637e-005 


1.948e-005 


3.941e-005 


2.115e-005 


128 


7.569e-006 


3.928e-006 


9.974e-006 


5.221e-006 


256 


1.340e-006 


1.340e-006 


2.482e-006 


1.292e-006 


512 


3.102e-007 


6.738e-007 


5.919e-007 


3.208e-007 



Table 1: u n+l and w™ +n are FEM solutions on T n+1 and Shishkin mesh T s n+n for Example 16. II 
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Let e = 1CT 10 . Table 121 shows the accuracy of Qi, the intersections of u s and u s+l . 
Denote x- and y- coordinates of Qi by x(Qi) and y{Qi), respectively. Note that Qi has 
better accuracy than u 8+1 . The reason is yet to be discovered; see Table El 



i 


x(Qi) 


\y(Qi)-«Qi))\ 


\y(Qi)-u*MQi))\ 


2 


0.2499999996000000 


7.499999579718697e-011 


4.999999719812465e-011 


3 


0.2500000004000000 


2.500008533523612e-011 


8.326672684688674e-017 


4 


0.4999999992000000 


3.500000012035542e-010 


2.999999970665357e-010 


5 


0.5000000008000000 


5.000011515932101e-011 


1.110223024625157e-016 


6 


0.7499999988000000 


6.625580639685325e-009 


6.700580590379701e-009 


7 


0.7500000012000000 


7.500006171667906e-011 


1.110223024625157e-016 



Table 2: errors at Qi with e — 10 10 on T 8 and T 8+1 for Example 16.11 

Plotted in FigureEJare the convergence curves in the maximum norm — M n+1 ||oo,[o,x n ] 
for e = 10~ 5 and e = 10~ 10 , respectively. They clearly indicate the convergence rate is 
proportional to rT 2 . It verifies Remark 14. 4| see Figure |2j 

Example 6.2. We examine the problem of a reaction diffusion equation as another example 
of©. 

-eu"(x) + u{x) = x, u(0) = u(l) = 0. (44) 

The exact solution is 

u(x) =x _ e _ 2/v - (45) 

The exact solution u(-) has boundary layer at x = 1, and is nearly linear outside the 
boundary layer. Also, reaction diffusion equation has relatively stable matrix A compared 
with convection diffusion equation. Due to these reasons, the FEM solutions of (|44|) is better 
than the FEM solutions of (EE2l . 

For e = 10~ 10 , we compute the FEM solution u 4 and u A+1 on the grid T 4 and T 4+1 , 
where T 4 is uniform mesh on [0, 1] and T 4+1 is modified by adding one point at the center 
of last interval; Figure H3-a. 

By adding one point si G (x n , 1) with Si — x n = y6e as in Lemma 14. 1[ we use new 
grid 7~ n+1 , and denote its FEM solution as m 4+1 . As shown in Figure HU-b, u 4+l is almost 
overlapped with uj, the interpolation of exact solution; see Figure HHb. 
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Let 9 = min{|, y ^^^}. We construct shishkin mesh T s n+n by dividing [0,1 — 9] and 
[1 — 9, 1] into n equidistant subintervals. Table El present the errors of u n+1 . Compared with 
m" + ", the FEM solutions using Shishkin mesh T™ +n , the errors are smaller and e-uniform. We 
omit the convergence curve and error table of Qi, since all those errors are within computer 
errors (around 10~ 14 ). 




e = 


: 10" 5 


e - 


= 10~ 10 


n 


II n+l f.n+l\\ 

||% —U ||oo,[0,x„] 


||% — U s ||oo,[0,Xn] 


II n+l f n+l || r 


|L,«+1 n,n+n\\ 

n] 11% -% l|oO,[0,X„] 


4 


1.665e-016 


1.517e-004 


1.110e-016 


4.980e-007 


8 


1.110e-016 


5.415e-005 


2.220e-016 


1.868e-007 


16 


2.220e-016 


2.161e-005 


3.331e-016 


8.451e-008 


32 


2.220e-016 


7.391e-006 


3.331e-016 


4.054e-008 


64 


3.331e-016 


1.300e-006 


5.551e-016 


1.984e-008 


128 


4.441e-016 


1.159e-009 


5.551e-016 


5.551e-016 


256 


2.459e-013 


2.948e-007 


6.661e-016 


9.795e-009 


512 


5.440e-015 


2.865e-007 


7.772e-016 


4.847e-009 



Table 3: u n+1 and u™ +n are the FEM solutions on T n+1 and Shishkin mesh T™ +n for Exam- 
ple E2 



Example 6.3. This example presents a demonstration of Corollary 13. 51 Using the FEM, we 
aim to find the Green function ( solution of) 

-eV + u = <? a , u(0) = 0,«(1) = 0, (46) 

where 5 a is delta function with peak at a G (0, 1). Denote a function as 

g {x) = — -. (47) 

° v ' po/e p— a/e v ' 



The exact solution of (jjBj) is 



K a g(x) < x < a 

U =\ Ka9 («£>) a< X <l (48) 



where K a ~ a/e is a constant depend on a and e. 
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T 6+1 used in Figure EJ-a is modified from T 6 by adding one point s\ at the center of the 
last interval, while T 6+1 used in Figure EJ-b is modified by adding si at the center of (x 3 , x 4 ). 
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(a) u 6 and u 6+1 (dotted lines) for f < a < 1. (b) u 6 and it 6+1 (dotted lines) for | < a < |. 



Figure 4: FEMs with e = 10~ 5 for Example IO 



< a < 1 



| < a < | 



2 .1714 

3 .3158 

4 .4588 

5 .6016 

6 .7445 



i x (Qi) \y{Qi) ~ u(x(Qi))\ x(Qi) \y(Qi) - u(x(Qi))\ 



1.7347e-018 
6.9389e-018 
2.7756e-017 







.1714 
.3158 

.6842 
.8286 



2.7756e-017 
1.1102e-016 

1.7764e-015 
4.4409e-016 



Table 4: errors at Qi for e = 10 5 for Example 16.31 



7 Further Remarks 

This paper is devoted to finite element methods for singularly perturbed boundary value 
problems. An interesting behavior is discovered: One can add arbitrary many points in 
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one of the grids, while the corresponding FEM solutions always have the common inter- 
sections {Qi} in all other intervals. Moreover, a practical and efficient ^-uniform mesh is 
developed. The FEM solution under this mesh can be viewed as a non-singularly perturbed 
BVP perturbation problem, and all general FEM error analysis can be applied. 

In both Example 16.21 and Example 16. 3[ the errors are within computer error. However, 
the errors of Example 16.11 is visible errors relative to computer error. The main reason is 
the exact solution of Example 16. II is almost quadratic, while our approximation is based on 
linear finite element space. To increase accuracy, one can generalize the results to the higher 
order finite element space. If the exact solution has several boundary layers, it can also be 
generalized to isolate each boundary layer. 

Although the exact solution of Example 16 .11 is nearly quadratic, the accuracy of intersec- 
tions {Qi} is almost within computer error. We know u n+1 has the accuracy of w n+1+0 °, while 
{Qi} has the accuracy of u n+0 °. The only difference of two is the interval (x n ,Si) of width 
0(y/e) or 0(y/e). In fact, this causes the error difference from Example 16.11 It might be 
interesting to discover the reason behind. It leads to the error analysis of non-quasiuniform 
meshes. 

It is very challenged to generalize the idea to isolate boundary layer in higher dimensional 
cases. On the other hand, Lemma f5. II provided a necessary and sufficient condition to verify 
the behavior of oscillation of specific FEM solution. However, it is not handy enough to 
explain why the oscillation behavior is common to FEM solutions. In general, the problem 
of determining in what cases the FEM solutions will or will not oscillate remains open. 
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